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g : ABSTRACT 

<N : 

^ , We present a detailed analysis of XMM-Newton EPIC-pn data for the Seyfert- 

^ '. 1 galaxy NGC 4593. We discuss the X-ray spectral properties of this source as 

! well as its variations with time. The 0.5 — 10 keV spectrum shows significant com- 

I— i! plexity beyond a simple power- law form, with clear evidence existing for a "soft 

p ^! excess" as well as absorption by highly ionized plasma (a warm absorber) within 

Q ! the central engine of this active galactic nucleus. We show that the soft excess is 

^ ! best described as originating from thermal Comptonization by plasma that is ap- 

d . preciably cooler than the primary X-ray emitting plasma; we find that the form 

! of the soft excess cannot be reproduced adequately by reflection from an ionized 

^ i accretion disk. The only measurable deviation from the power-law continuum in 

! the hard spectrum comes from the presence of cold and ionized fluorescent iron 

00 

Q I Ka emission lines at 6.4 and 6.97 keV, respectively. While constraints on the 

^ ! ionized iron line are weak, the cold line is found to be narrow at CCD-resolution 

^ I with a flux that does not track the temporal changes in the underlying contin- 

I uum, implying an origin in the outer radii of the accretion disk or the putative 

^ I molecular torus of Seyfert unification schemes. The X-ray continuum itself varies 

. !_i I on all accessible time scales. We detect a ~ 230 s time-lag between soft and hard 

^ I EPIC-pn bands that, if interpreted as scattering timescales within a Comptoniz- 

5^ \ ing disk corona, can be used to constrain the physical size of the primary X-ray 

source to a characteristic length scale of ~ 2rg. Taken together, the small implied 
coronal size and the large implied iron line emitting region indicate a departure 
from the current picture of a "typical" AGN geometry. 
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1. Introduction 

It is tliought tliat supermassive black holes are the central engines of active galactic 
nuclei, and may inhabit the cores of most galaxies (Salpeter 1964, Zel'Dovich 1965, Lynden- 
Bell 1969). Accretion disks around such objects provide some of the best natural laboratories 
for examining the effects of an extreme gravity environment on infalling material. With the 
enhanced throughput and spectral resolution of modern X-ray telescopes such as Chandra, 
XMM-Newton and Suzaku, it has now become possible to observe these effects with the 
precision necessary to make robust measurements of the innermost properties and overall 
structures of accretion disks. Evidence for relativistic emission line broadening, spectral and 
temporal variability over a wide range of time scales, and photoionized (warm) absorption 
from outflowing highly-ionized gas along the line of sight are only some of the interesting 
spectral detections being made in many active galaxies. 

In this paper, we present an X-ray spectral analysis of the central regions of the canonical 
Seyfert-1 galaxy, NGC 4593. Galaxies of this type are of particular interest for the study of 
accretion onto supermassive black holes since they are typically oriented such that we can 
view the accreting black hole free of substantial obscuration or absorption from surrounding 
circumnuclear material. Furthermore, we believe that significant amounts of the continuum 
radiation seen in radio-quiet Sy-1 galaxies originate in the disk rather than, for example, 
a relativistically beamed jet (as is the case for BL-Lac objects). As such, we are better 
positioned to study the inner parts of the accretion disk in Sy-1 galaxies than in other 
systems. 

NGC 4593 is a spiral galaxy with a central bar classified as Hubble type SBb. At 
a redshift of z = 0.009, it lies at a proper distance of 37.9 Mpc toward the constellation 
Virgo. This corresponds to an angular size distance of 37.6 Mpc, and a luminosity distance 
of 38.3 Mpc using Hq = 71kms~^ Mpc~^, fim = 0.27, and f^A = 0.730 The galaxy has an 
apparent visual magnitude of 11.67 and an approximate angular diameter of 3.9 x 2.9 arcmin. 
As already noted, it hosts an active galactic nucleus (AGN) of a Sy-1 type (Lewis, MacAlpine 
& Koski 1978). Previous studies of the source with EXOSAT demonstrated a soft excess 
(Pounds & Turner 1988), and B eppoS AX daia, confirm a broad absorption dip of 15% below 
1 keV which may be attributable to the presence of a warm absorber along the line of sight 
(Kaastra & Steenbrugge 2001). ASCA spectra display a slightly broadened cold iron line at 
6.4 keV, in addition to evidence for a warm absorber within the system (Nandra et al. 1997; 
Reynolds 1997). The source also displays significant variability in fiux. Between two ASCA 
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observations 3.5 years apart, the 2 — lOkeV flux of tliis source increased by ~ 25%, tliougli 
no significant variability of the iron line was detected between the two pointings (Weaver et 
al. 2001). Within the ASCA observation, Reynolds (1997) witnessed a count rate decrease of 
~ 60% in 10 ks, with smaller flares and dips throughout the data set. In terms of the overall 
properties of the system, Reynolds (1997) calculated a luminosity of Lx — 8.53 x 10^^ ergs~^ 
(2 — lOkeV) with ASCA, and more recently, McKernan et al. (2003) used Chandra to derive 
a luminosity Lx = 5.37 x lO^^ergs"''' (2 — lOkcV). It should be noted that this change in 
Lx is larger than the flux calibration uncertainty of the observation, so this does appear to 
be a robust flnding. 

In this paper, we present a 76 ks exposure of NGC 4593 with the X MM- Newton/ EPIC- 
pn instrument from 2002 June 23/24. We discuss our observations of the time-averaged 
2 — lOkcV spectrum and cold iron line in more detail in Reynolds et al. (2004; hereafter 
R04). Here, we will examine the spectral and temporal variability of NGC 4593, as well as 
the soft X-ray spectrum. §2 discusses the method and software employed for the reduction 
of the data referred to in this work. Spectral fitting analysis is detailed for the EPIC-pn 
instrument in §3, and variability studies are examined in §4. We present our discussion in 
§5, including a comparison of our results with those pubhshed in earlier studies of NGC 4593 
and a discussion of the disparate spatial scales implied for the X-ray emitting region and the 
iron fluorescence region. We also draw our conclusions in §5. 



2. Data Reduction 

In this paper, we use data taken with the European Photon Imaging Camera pn (EPIC- 
pn) camera. The data were obtained during revolution 465 of XMM-Newton, during which 
the pn was operated in its small-window mode to prevent photon pile up, using the medium 
fllter to avoid optical light contamination. The EPIC MOS-1 camera took data in the fast 
uncompressed timing mode, and the MOS-2 camera operated in prime partial W2 imaging 
mode. Although the MOS results will not be discussed further here, they mirrored the 
EPIC-pn data within the expected errors of calibration effects. The average EPIC-pn count 
rate for this source was 29.78 ctss~^. 

Data were also collected during the observation by the Reflection Grating Spectrometer 
(RGS). These data were examined but were deemed to be unhelpful for the analysis presented 
in this paper. Firstly, the RGS spectra have insufficient signal-to-noisc to allow a detailed 
study of the soft X-ray emission and absorption lines from the warm absorber in this source. 
They also contain insufficient counts to allow variability analysis. Finally, our examination 
of the RGS spectra have demonstrated that problems with the broad-band effective area 
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calibration of the RGS (which is strongly affected by a build-up of carbon on the CCDs) 
compromises any RGS-based results on the shape of the soft excess. For a detailed discussion 
of these calibration uncertainties, see Kirsch et al. (2005)@ For these reasons, will shall not 
consider RGS data further. 

The pipeline data were reprocessed using the Science Analysis Software and the cor- 
responding calibration files, version 6.5.0. From these, we rebuilt the calibration index file 
using cif build. For the EPIC-pn data the event files were mildly edited in spectral coverage 
to observe the region from 0.2 — 15keV, and bad pixels and cosmic ray spikes were removed 
via narrow time filtering using the evselect task within the SAS. No background flares were 
detected during the observation. Extraction of spectra followed the procedure used by Wilms 
et al. (2001), in which source and background spectra were generated using the xmmselect 
task. Response matrices and ancillary response files were created using rmf gen and arf gen, 
and the data were then grouped using the grppha task with a binning factor of 25ctsbin^^. 
Binning is required in order to get sufficient counts per bin to make spectral fitting a 
valid statistical process. Spectral modeling and analysis from 0.5 — lOkeV was performed 
using the XSPEC package version 11.3.2. (Arnaud 1996). Timing studies were performed 
using various routines in the XRONOS package (Stella & Angelini 1991). 

We have used the SAS epatplot task to compute the fraction of single, double, triple 
and quadruple events as a function of energy and compared these fractions to their nominal 
values as measured from weak source observations. For sources that are affected by pile up, 
these fractions deviate from the nominal values due to the higher probability of wrong pattern 
classification. No significant deviation from the nominal single and double distributions was 
found, indicating that our EPIC-pn observation of NGC 4593 is not affected by pile up. 



3. Spectral Analysis 

As shown in R04, the 2-10 keV EPIC-pn spectrum of NGC 4593 is very well described 
by a power-law with the exception of two emission features identified as the fluorescent 
Ka emission line of cold iron at 6.4 keV and the Lya recombination line of hydrogen-like 
iron at 6.97 keV (Fig. 1). The cold iron line likely arises from fluorescence on the surface 
layers of the outer accretion disk, optically-thick optical broad emission line clouds, or the 
putative "molecular torus" in response to hard X-ray irradiation from a hot corona in the 
inner accretion disk (Basko 1978; Guilbert & Rees 1988; Lightman & White 1988; George 



^See also http: / /xmm.vilspa.esa.es/external/xmin^w_cal/calib/index.shtml for a thorough discussion of 
these cross-cahbration issues. 



- 5 - 



& Fabian 1991; Matt et al. 1991). The ionized line, by contrast, could be formed either by 
ionized disk irradiation or by radiative recombination in highly-ionized outflowing material 
above the plane of the disk. 

However, extrapolating below 2keV reveals a clear soft excess as well as the effects of 
the warm absorber (Fig. 2; the "notch" in the range 0.7-1.0 keV signals the presence of the 
warm absorber). Since the basic origin of such soft excesses are still the subject of much 
debate, it is interesting to characterize the soft excess in as much detail as possible. This 
necessitates, however, a detailed modeling of the warm absorber. 

We model the effect of the warm absorber using the XSTAR photoionization code (version 
2.1kn3; originally developed by Kallman & Krolik 1995). In detail, we construct grids of 
models describing absorption by photoionized slabs of matter with column density and 
ionization parameter ^, 



where Li is the luminosity above the hydrogen Lyman limit, Ue is the electron number 
density of the plasma and r is the distance from the (point) source is ionizing luminosity. 
We have constructed 20 x 20 grids of models uniformly sampling the (log A^h, log^) plane in 
the range A^h = 10^° lO^^cm"^ and ^ = 1 — > lO^ergcm"^ s~^ . While these were made to 
be multiplicative absorption models and hence can be applied to any emission spectrum, the 
ionization balance was solved assuming a power-law ionizing spectrum with a photon index 
of r = 2, cut off at an energy Ec — 20keV. This is a good approximation of the typical 
AGN continuum. 

As a first, purely phenomenalogical attempt to describe the soft excess component, we 
have employed a single temperature thermal bremsstrahlung emission model (zbremss). To 
adequately describe the spectrum, we were required to include two warm absorbing zones 
as well as an iron-Ls edge from iron atoms in line-of-sight dust. A natural interpretation 
is that the WA has distinct components at physically different distances from the central 
engine. We note that attempting to add a third zone into the model results in no statistical 
improvement in fit. Including this absorption structure provides a good description of the 
data (x^ = 1813/1745; Fig. 3 and Table 1). We shall refer to this model as Model-1. The best 
fit value for the temperature characterizing this soft excess component is kT 0.21^q keV 
and it has a 0.5-10 keV flux of 2.011°;?? x lO-^^ergcm'^ s^^ . 

Of course, the description of the soft excess in Model-1 is purely phenomenalogical in 
nature. Ideally, we would like to describe the soft excess with a model that has a direct 
physical interpretation. Initially, we attempt to describe the soft excess as soft X-ray reflec- 
tion from a mildly-ionized accretion disk (see Brenneman & Reynolds 2006 for a successful 
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Fig. 1. — The 2 — lOkeV spectrum of NGC 4593 fit with a simple photoabsorbed power-law 
(phabs po). Note the residual iron features at 6.4 and 6.97keV. xV^of = 1879/1450(1.30). 
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Fig. 2. — The 0.5 — 10.0 keV spectrum of NGC 4593 fit with the photoabsorbed power-law 
from Fig. 1 and two Gaussians to model the iron lines. Note the clear evidence for a soft 
excess, possibly complicated by absorption features from a "warm absorber" within NGC 
4593. For this fit, xVdof = 11423/1736(6.58). 
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Fig. 3. — (a) The best fit for Model 1, including a zbremss component to represent the soft 
emission: x^/dof = 1813/1745(1.04). (b) The relative strength of the model components 
for Model 1. The 6.4 and 6.97 keV Gaussians are shown in red and green, respectively. The 
zbremss soft emission is in dark blue, and the photoabsorbed power-law is in light blue. 
Other absorption components are indicated in black. 
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application of this model to the Seyfert galaxy MCG-6-30-15). Operationally, we replace the 
thermal bremsstrahlung component in Model 1 with the ionized disk model ref lion created 
by Ross & Fabian (2005). Because the irradiated matter is also responsible for producing 
the Fe-Ka line in many sources, this model has the potential advantage of self-consistently 
describing the soft excess as well as the 6.4 keV emission line feature. 

Interestingly, we find that we cannot successfully describe both the soft emission and the 
Fe-Ka line with the disk reflection model. The resulting best fit is xV^of = 2652/1745(1.52), 
and we found that significant residuals remained on the soft and hard ends of the spectra. 
Even adding the Gaussian component back to exphcitly model the 6.4 keV iron line (and 
hence removing the constraint that the disk refiection must describe this feature), the disk 
reflection model was unable to adequately fit the form of the soft excess. We thus conclude 
that ionized disk reflection is not an important process in shaping the soft X-ray spectrum 
of NGC 4593. Given the relative narrowness of the Fe-Ka line in this source we believe it is 
not originating in the inner disk, and hence the above is not a surprising conclusion. 

Alternatively, we can postulate that the soft excess is due to thermal Comptonization 
by plasma at a temperature between that of the disk and the hard X-ray corona (possibly in 
a transition zone). Operationally, we use the comptt model (Titarchuk 1994) to model the 
soft excess. Hereafter, we shall refer to this spectral model as Model 2. This preserves the 
physical realism of the soft emission arising from the accretion disk, but allows the iron line 
to be produced elsewhere, perhaps farther away in the molecular torus where it would not be 
as broad by nature. Model 2 reaches a best fit of xV^of = 1759/1744(1.01). The best fit and 
model for Model 2 is shown in Fig. 4. Here we have frozen the seed photon temperature at 
Tq — 50 eV and have kept a slab geometry for the corona. The best fit Comptonizing plasma 
temperature is kT^ = 42keV with an optical depth of r = 0.12, though both parameters 
were not very well constrained by the fit. This is not surprising, given that both arc equally 
involved in shaping the spectrum via the Compton-y parameter: y oc rTc. We note that 
for Model 2 the equivalent widths of the 6.4 and 6.97 keV Gaussians remain approximately 
unchanged from their values in the Model 1 fit. 

Table 1 reports the parameter values and goodness of fit for both the final best fit- 
ting phenomenalogical model (Model 1, with zbremss and including all other components 
added in), and the more physical model (Model 2, including comptt and all other compo- 
nents). Using the Model 2 fit from Table 1, the total 0.5 - lOkeV fiux is Fx = 6.74 x 
10~^^ergcm~^s~^ . Assuming a fiat universe WMAP cosmology, this corresponds to a lu- 
minosity of Lx = 1.21 X lO^^ergs"^ . Considering only the energy range from 2 — lOkeV, 
-^2-10 = 7.40 X lO^^ergs"^. This is roughly 21% greater than the 2 — lOkeV luminos- 
ity observed by McKernan et al. (2003), but only about 86% of the value from the ASCA 
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Fig. 4. — (a) The best fit for Model 2, including a comptt component to represent tlie soft 
emission in place of the zbremss component of Model 1. x^/dof = 1808/1744. (b) The 
relative strength of the model components for Model 2. The color scheme is the same as in 
Fig. 3b. 
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observation reported by Reynolds (1997). 

The warm absorber modeled with oTir two created xstar tables suggests a multi-zone 
structure in column density and ionization, as discussed above. Though our model does not 
parameterize the covering fraction of the absorbing gas, we can approximate this fraction 
by representing the warm absorber with the Ovil and OviII edges at 0.74 and 0.87kcV, 
respectively, as has been done by many other authors examining warm absorbers (e.g., 
Reynolds 1997; McKernan et al. 2003). In so doing, we find that to 90% confidence, the O Vll 
edge has an optical depth of r = O.IGIqIqi whereas the Ovill edge is much less prominent 
at T < 0.01. We therefore chose to use the Ovil edge as a diagnostic for determining 
the covering fraction of the absorbing gas in NGC 4593. If the Ovil edge accounts for 
the absorption by this gas, the photons absorbed should be reprocessed and radiated in 
emission as the O VII radiative recombination line at 0.574 keV. Assuming that this line will 
be narrow, we fit a Gaussian to the data. Adding in this component bettered the overall 
goodness-of-fit by A^^/Adof = —53/ — 1 as compared to the two-cdgc model for the warm 
absorber, showing that this line is statistically robust. The Ovii RRC line had an upper 
hmit equivalent width of EW < 6.35 eV and a flux of 2.82 x 10"^ phcm~^ s~^ , such that 
the line contributed ~ 7182 total photons to the model over the total observation time. 
The Ovil edge, by contrast, contributed ~ 31991 total photons. Taking the ratio of the 
reprocessed photons to the absorbed photons, we can estimate a covering fraction for the 
warm absorber: / « 0.23. 

4. Variable Nature of the Source 

NGC 4593 displays significant variability (over a factor of two in 0.5 — lOkeV flux) 
during the course of this observation. In this Section, we examine the detailed variability 
properties of this source. 

4.1. Spectral Variability 

We have analyzed the spectral variabihty of the source as seen in the EPIC-pn data. As 
an initial assessment of spectral variability. Fig. 5 plots the data-to-model ratios of consecu- 
tive 10 ks segments of the EPIC-pn data using the comptt best flt spectrum model discussed 

in §3. Experimentation suggested that intervals smaller than 10 ks contain insufficient counts 
to maintain the integrity of the spectrum. As well as the obvious variations in the normal- 
ization of the spectrum, there are clear changes in slope between segments (with the source 
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Table 1: Best fit parameters for the EPIC-pn spectrum 
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A riA+1 46 1 , 1 a22 

9.29_g;^5 X 10^^ 

9 7P;+0:i2 
^- ' '-'-0.32 


WA 2 


A^H2 ( cm-2 ) 
log ^2 (erg cm" 1 
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1 Ol +0 02 


1 Ar>+0 04 

16.921^:^1 


po 


r 




1 87 -1- n nn/i 
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'^-0.03 




flux ( erg cm"^ s" 




5.11^°;°^ X 10-13 


4.44l°:^i X 10-13 


zgauss 


E(keV) 




6.4 


6.4 




a(keV) 




0.10 ±0.08 


0.10 ±0.01 




flux (ergcm""^s" 




4.46 ± 0.32 X 10-15 


4.i5i^:i X 10-^' 




E{keV) 




6.97 


6.97 




a(keV) 




0.10 ±0.02 


0.10 ±0.06 




flux ( erg cm"^ s" 




9.47^1°? X 10-16 


7.7711^^ X 10-16 


zbremss 


A;T(keV) 




n 91+0.003 
'-'•^J^-0.005 






flux ( erg cm"^ s" 




2.01^°i3 X 10-1^ 




comptt 


To(keV) 
A;T(keV) 






0.05 

1 2+°°6 




flux (ergcm^^s" 






6.46l°|° X 10-1^ 


xVdof 




1813/1745(1.04) 


1759/1744(1.01) 



The energy range from 0.5 — 10.0 keV is considered for the EPIC-pn. Best fit Model 1 contains a 
zbremss component to parameterize the soft excess below ~ 2 keV, while best fit Model 2 represents 
this component with a comptt model. All quoted error bars are at the 90% confidence level. All 
redshifts used were frozen at the cosmological value for NGC 4593: z = 0.009. 
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becoming softer as it brightens). There is also a feature of variable equivalent width seen at 
~ 6.4 keV, corresponding to the cold Fe-Ka line. Interestingly, though, we see no variable 
discrete features in the soft X-ray spectrum suggesting that there is no significant warm 
absorber variability during our observations. 

To examine variability of the iron lines in more detail, Fig. 6 plots the renormalized 
interval data over the time-averaged data from 2 — 10 keV. Direct spectral fitting of the the 
cold and ionized Fe-K hnes in each interval with Gaussian profiles illuminates the nature of 
their variabihty (or lack thereof). As can be seen in Table 2 and Figs. 7, the cold line exhibits 
an approximately constant flux and a variable equivalent width. The constant equivalent 
width hypothesis is rejected with greater than 99% significance, whereas a constant flux 
hypothesis is consistent with the data at the 20% level. Thus the evidence suggests that the 
cold iron line flux does not appear to respond to changes in the X-ray continuum. Constraints 
on the ionized line are not strong enough to rule out either a constant flux and equivalent 
width, as is shown in Fig. 8 and Table 2. 

The simplest way to interpret the lack of continuum response of the cold iron line is 
to suppose that it exists on spatial scales with light travel times greater than the duration 
of our observation. For an observation time of 76 ks, as is the case here, the light would 
travel approximately 150 AU, or roughly 1.9 x 10^ if NGC 4593 harbors a 8.1 x 10^ Mq 
black hole at its core (see discussion at the end of §4.2). The lack of response for the iron 
line flux would suggest that this size is a lower limit for the size of the line emitting source. 
This would place the line emission in the outer regions of the accretion disk or the putative 
molecular torus of the Seyfert unification scheme. 



4.2. Continuum Power Spectrum and Spectral Lags 

The hght curve for the full-band (0.5-10.0 keV) data set is shown in Fig. 9a. The hght 
curve is characterized by an initial rapid drop up to ~ 5 ks, followed by a relatively steady 
increase peaking at ~ 45 ks, a drop until ~ 65 ks, and then a final increase that shows possible 
signs of tapering off around the end of the observation at 76 ks. There are indications of 
smaller scale variability on time scales of < 1000 s, with slight flares and drops occurring 
throughout the data set. 

Fig. 9b shows the Leahy normahzed power spectrum of the 0.5 — lOkeV EPIC-pn light 
curve. The power spectrum demonstrates that the variability of the source becomes domi- 
nated by Poisson statistics on timescales of about 200 s once it drops to a power of ~ 2, as 
expected. Before this point, the slope of the power spectrum for low frequencies is — 2.63lo;25- 
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Fig. 6. — The cold and ionized Fe-Ka lines from the time-separated EPIC-pn spectrum. 
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Table 2: The cold and ionized iron lines in the EPlC-pn spectrum of NGC 4593 



Line 


FWHM ( km s^^ ) 


Flux ( erg cm ^ s ^ ) 


EW ( eV) 


Cold Fe-Ka 
(6.40 keV) 


13245 ± 2208 
8831 ± 2208 
9934 ± 2208 

13246 ±3311 

8831 ± 2208 
11038 ±4415 
8831 ± 2208 
11038 ±3311 


7.41 ± 0.86 X 10-13 
4.98 ± 0.89 X 10-13 
5.76 ± 0.85 X 10-13 

7.26 ±1.04 X 10-13 

5.27 ±0.86 X 10-13 
4.08 ±0.93 X 10-13 
5.61 ±0.88 X 10-13 
5.82 ± 1.17 X 10-13 


232 ± 27 
126 ± 22 
134 ± 19.8 
153 ±21.8 
105 ± 17.2 
88.5 ±20.2 
137 ±21.5 
124 ±24.9 


Ionized Fe-Ka 
(6.97 keV) 


7095 ± 6081 
50678 ±83111 
6081 ± 6081 
14190 ± 11149 
9122 ± 6081 
18244 ± 9122 
24325 ± 9122 


1.97 ± 1.41 X 10-13 

2.89 ± 5.05 X 10-13 
2.09 ±1.54 X 10-13 
2.04 ± 2.22 X 10-13 
2.52 ± 1.85 X 10-13 
3.06 ±3.48 X 10-13 

4.90 ± 3.00 X 10-13 


64.3 ±46.0 
80.3 ± 140 
51.9 ±38.3 

45.6 ±49.6 

53.7 ±39.4 
71.3 ±57.8 
132 ± 80.9 









Each row represents a time interval of 10 ks in the observation. Fits were done in XSPEC 
using a phabs po model for the continuum and two Gaussian lines with parameters fit to 
the data in each time interval. Error bars are at la. The blank lines in the ionized Fe-Ka 
table represent time intervals for which a robust fit to the data could not be achieved. 
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Fig. 7. — Variation of the flux (a) and equivalent width (b) of the cold Fe-Ka line (6.4 keV) in 
10 ks intervals over the course of the 76 ks observation. A fit to the data with a constant flux 
model yields x^/dof =10/7 (1.43) implying a 20% probability that the data are consistent 
with the constant flux model. The constant equivalent width fit yields x^/dof = 22/7 (3.14) 
allowing the constant equivalent width hypothesis to be rejected at more than the 99% 
significance level. 
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Fig. 8. — Variation of the flux (a) and equivalent width (b) of the ionized Fe-Ka hne 
(6.97 keV) over the course of the observation. Here the uncertainty in the data is so great 
that a constant flux cannot be ruled out: x^/dof = 3/7 (0.43). A constant equivalent width 
is similarly likely: x^/dof = 3/7 (0.43). 
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This value for the PSD slope is consistent with results from several other Sy-1 AGN such 
as NGC 3783 (Markowitz 2005), MCG-6-30-15 (Papadakis et al. 2005; Vaughan et al. 2003), 
NGC 4051 (at high frequencies; McHardy et al. 2004), and Mrk 766 (also at high frequencies; 
Vaughan & Fabian 2003). 

Given the range of X-ray energies available in this dataset, it is also useful to study 
the continuum variability in different energy bands. Hardness ratios and time lags between 
different energy bands in the continuum, in particular, can help constrain the emission 
mechanisms of the source as well as the physical scale of the emitting region(s). We divide 
the X-ray spectrum into three energy bands, each with an approximately equal number of 
counts (~ 5 X 10^). The soft band ranges from 0.5 — 0.84 keV, the mid band extends from 
0.84 — 1.5 keV, and the hard band covers 1.5 — 10.0 keV. The light curves for all three bands 
are represented in Fig. 10. In order to test whether a time lag exists between these energy 
bands, we compute the cross-correlation function between the three energy bands using the 
XRONOS package (Stella & Angelini 1991) as well as the discrete correlation function (DCF; 
Edelson & Krolik 1988). We measured the hard-to-soft band lag to be 230 ± 50s (i.e., the 
hard band lags the soft band by this amount of time; Fig. 11). 

The most straightforward interpretation of these time delays envisions them as repre- 
senting the finite scattering time within the Comptonizing disk-corona system that is thought 
to be responsible for the primary X-ray production. One should be cautioned, however, that 
this interpretation invokes numerous underlying simplifications with regard to the physical 
mechanism that produces these types of time lags between energy bands. See Pottschmidt et 
al. (2003) and Nowak et al. (1999, 2002) for more detailed accounts of other factors that can 
be responsible for or affect the observed time lag between the hard and soft energy bands. 
We assume here that both the hard and soft X-ray photons originate from seed UV photons 
from the disk, which are then Comptonized by a corona of relativistic electrons surrounding 
the disk. The energy of a given photon in our simplified scheme will thus depend on the 
number of scatterings it undergoes with the electrons, and one may therefore infer that the 
more energetic photons have experienced more interactions. It is thus possible to estimate 
the size of the coronal region by measuring the time lag between the peaks of the hard and 
soft photon light curves of the source, as we have done, and making the simplified assumption 
that the corona possesses a slab or spherical geometry. With each scattering, a seed photon 
gains on average a fractional energy is.E / E ^ -iksT /uieC^ , where T is the temperature of 
the electrons in the corona. So the energy of a Comptonized photon after a given number 
of scatterings will be: 



where Eq is the initial seed photon energy and n is the number of scatterings (Rybicki & 




(2) 
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a) 




Frequency (Hz) 



Fig. 9. — The light curve (a) and power spectrum (b) for the time-averaged pn data from 
NGC 4593. Variations on timescales as small as hundreds of seconds appear visible in the 
source. Note that the frequency at which the power spectrum flattens into Poisson noise is 
about 5 X 10"'^ Hz. Inverting this means that the smallest timcscalc of variability wc can 
reliably observe from this source is ~ 200 s. At low frequencies, the slope of the power 
spectral density curve is — 2.631;q;25. 
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Observation Time (s) 



Fig. 10. — Light curves for tlie different energy bands: the soft band runs from 0.5 — 0.84 keV 
(solid black), the mid band runs from 0.84 — 1.5 keV (dashed red), and the hard band runs 
from 1.5 — 10.0 keV (dotted green). 
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Lightman 1979; Chiang et al. 2000). The time delay of this photon relative to the seed 
photon source will be roughly proportional to the number of scatterings, t ~ ntQ. Here 
Iq ~ It/c where It is the mean free path for Thomson scattering in an optically-thick corona 
or, if the corona is optically-thin, simply the size of the corona. We take the effective photon 
energy in the soft range (0.5 — 0.84 keV) to be 0.67 keV. For the hard range (1.5 — lOkeV) 
we calculate the effective photon energy to be 2.17 keV. Both estimates are based on taking 
weighted averages (based on flux) of the energies in the soft and hard ranges, and accounting 
for the energy-dependent effective area of the pn. Assuming that the energy of the coronal 
electrons ranges from 50 — 100 keV, our measured soft-hard time lag of ~ 226 s suggests that 
the corona occupies a region around the disk 1.91 — 3.34 x 10^^ cm in size. This is derived 
assuming the size of the corona is given by 



We will assume that the central black hole in NGC 4593 has a mass of 8.1 x 10^ Mq based 
on reverberation mapping (Nelson & Whittle 1995; Gebhardt et al. 2000). Using this mass 
for the NGC 4593 black hole, the characteristic length scale of the corona would be in the 
range of 1.6 — 2.8rg. 



We have obtained a 76 ks XMM-Newton observation of the Sy-1 galaxy NGC 4593. 
An examination of the best fitting EPIC-pn spectrum shows that the continuum is well 
modeled by a photoabsorbed power-law with Nn = 1.97 x 10^° cm~^ , T = 1.75, and a flux of 
4.44 X 10~^^ erg cm~^ . The best fit to the hard spectrum can be achieved by including two 
Gaussian emission lines representing cold and ionized fluorescent iron at 6.4 and 6.97 keV, 
respectively. 

We see clear evidence for a complex warm absorber in the EPIC-pn spectrum. Fitting 
a grid of photoionization models computed using the XSTAR code, we infer a warm absorber 
with two physically and kinematically distinct zones: one with a column density of A^h = 
9.29 X lO^^cm"^ and an ionization parameter of log^ = 2.75, the other with Ah = 1-13 x 

10^^ cm~^ and log,^ = 1.70, which is likely more distant from the central engine. We robustly 
detect the L3 edge of neutral iron presumed to exist in the form of dust grains along the linc- 
of-sight to the central engine, as cited by McKernan et al. (2003). This edge has a column 
density of A^pe = 8.32 x lO^^cm'^ . 




(3) 



5. 



Discussion and Conclusions 



5.1. Summciry of Results 
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Fig. 11. — The cross-correlation function for the soft-to-hard energy bands in NGC 4593 
in red (with data points), plotted with the best fit curve in sohd black. The vertical line 
represents a zero time delay between the two bands; the oflF-centeredness of the curve peak 
indicates that the hard band lags the soft band by 226 ± 53 s at 90% confidence. 
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A soft excess below ~ 2keV is also seen in the data, and can be accounted for either 
phenomenalogically by a redshifted bremsstrahlung component (Model 1) or, more effectively 
and physically, a component of Comptonized emission from an accretion disk of seed photons 
upscattered by a plasma of relativistic electrons existing in either a corona or the base of a 
jet near the disk in some geometry (Model 2). The latter scenario yields the best statistical 
fit to the data, with a seed photon temperature of Tq = 50 eV, an electron temperature of 
kT = 42keV, a plasma optical depth of Tp = 0.12 and a flux of 6.46 x 10~^^ergcm~^s"^ . 

Lastly, cold and ionized Fe-K emission lines were wanted by the fit at 6.4 and 6.97 keV, 
respectively. The cold line Gaussian had a time-averaged equivalent width of EW — 131 eV, 
while the ionized component had an EW — 45 eV. Although the continuum varied with 
time over the course of our observation, the flux of the cold Fe-Ka showed marginal evidence 
for variability between successive 10 ks intervals. The equivalent width of this line, on the 
other hand, is shown to vary significantly. The simplest interpretation of this result is to 
suppose that the cold line originates from a region with a light-crossing time larger that the 
length of our observation. For a black hole mass of 8.1 x 10® Mq, this places the cold line 
emitting region beyond about 2000rg from the black hole, i.e., in the outer accretion disk 
or the putative molecular torus of Seyfert unification schemes. Our statistics on the ionized 
Fe-K line are insufficient to constrain its variability properties — these data are consistent 
with both constant fiux and constant equivalent width. 

We have detected a 226 s time-lag between the hard and soft EPIC-pn bands, with 
the hard band lagging the soft. In a simple model in which this corresponds to scattering 
times within a Comptonizing corona, we conclude that the corona can only possess a size of 
~ 1.6 - 2.8r^. 



5.2. Implications for the X-ray Emission Region 

The narrowness of the fiuorescent iron line together with the lack of response of this fine 
to changes in the hard X-ray continuum suggest an absence of a cold, optically-thick matter 
within the central ~ lO^r^ of the accretion disk. The standard framework for accommodating 

such a result is to postulate that the inner regions (r < lO^rg) of the accretion flow have 
entered into a radiatively-inefficient mode which is extremely hot (electron temperatures of 
T ~ 10^*^ K), optically-thin, and geometrically-thick (e.g., Rees 1982; Narayan & Yi 1994). 
Within this framework, the hard X-ray source is identified as thermal bremsstrahlung or 
Comptonization from this structure. However, the variability of the X-ray continuum is 
inconsistent with X-ray emission from a structure which is ~ lOOOrg in extent. The rapid 
continuum variability and the short time lag between the hard and soft X-ray photons dictate 
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that, at any given time, the X-ray emission region is only ^ 1 — Svg in extent. 

We are therefore led to consider alternative geometries and/or origins for the X-ray 
source. We consider three possibilities. Firstly, the inner accretion flow may indeed be 
radiatively-inefficicnt, but the observed X-ray emission may come from the compact base 
of a relativistic jet powered by black hole spin. Secondly, the X-ray emission may indeed 
originate from the body of the radiately-inefficient flow but, at any given instant in time, be 
dominated by very compact emission regions within the flow. Such emission regions may be 
arise naturally from the turbulent flow or, instead, may be related to a magnetic interaction 
between the inner parts of the flow and the central spinning black hole (Wilms et al. 2001; 
Ye et al. 2007). In this case the iron line could still plausibly originate at a distance of 
thousands of and thus be quite narrow. Finally, it is possible that the accretion disk is 
radiatively- efficient (and hence optically-thick and geometrically-thin) close to the black hole 
and supports a compact accretion disk corona. Of course, an immediate objection to this 
scenario is the lack of a broadened iron line. However, as demonstrated in R04, it is possible 
the iron line is so broad that it is buried in the noise of the continuum. Alternatively, the 
accretion disk surface may have an ionization state such that iron line photons are effectively 
trapped by resonant scattering and destroyed by the Auger effect. Longer observations with 
better spectral and timing resolution will be necessary in order to conflrm our results and 
differentiate between these different possible scenarios. 



5.3. Conclusions 

We have shown that the Sy-1 galaxy NGC 4593 has a continuum spectrum that is flt 
remarkably well by a simple photoabsorbed power-law above ~ 2 keV. Below this energy, 
we see evidence for spectral complexity that can be attributed to the presence of a possible 
multi-zone layer of absorbing material intrinsic to the source, as well as a soft excess that 
cannot be explained by a reflection model from an ionized disk. Also arguing against a disk- 
reflection-dominated source is that unlike other sources of its kind (e.g., MCG-6-30-15), 
NGC 4593 has relatively narrow cold and ionized Fe-Ka line features at 6.4 and 6.97 keV, 
respectively. See R04 for a more thorough discussion of these lines and their possible physical 
interpretations. We can say that the cold line is most likely formed either quite far out in 
the accretion disk, or possibly in the putative dusty torus region surrounding the central 
engine, based on their narrowness and lack of signiflcant variability over the 76 ks duration 
of our observation. We flnd no evidence for reflection features from the inner accretion disk 
(e.g., a "broad iron line") in the spectrum. 

Based on the time lag between the soft and hard spectral bands, we estimate that the 
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corona occupies the a region around the central source on the order of ~ 1.6 — 2.8rg, assuming 
that NGC 4593 harbors a 8.1 x 10^ black hole at its core and that the energies of the 
electrons in the corona range from ~ 50 — 100 keV. This estimate for the coronal size is 
reinforced by our measurements of continuum variability on timescales as small as ~ 200 s, 
equivalent to a light travel time of ~ 4rg for the black hole mass in question. 

Taken together, the implications of a narrow iron line emitted far out in the disk or 
torus and the small coronal size in NGC 4593 present an atypical picture of an AGN. We 
postulate that the primary X-ray source is associated with the compact base of a jet or 
compact emission regions within a much larger optically-thin accretion flow. Alternatively, 
the accretion disk may be radiatively-efficient with a compact corona but not display a broad 
iron line due to the effects of extreme broadening or disk ionization. 
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